Galaxy formation and evolution in the Millennium Simulation

This Thesis addresses the topic of galaxy formation and evolution in the universe. In collaboration with D. Croton, G. de Lucia, V. Springel, and S.D.M. White, I made use of the Millennium simulation, a very large N-body simulation of dark-matter evolution in a cosmological volume carried out at the MPA in 2005 by Springel 2005, to explore the predictions made by the most recent generation of semi-analytic models for galaxy formation. These models are incorporating a new mode of feedback from active galactic nuclei (AGN), which have their origins in super-massive black holes accreting mass and turning it into energy. Because of its observational signature in the radio regime this feedback is called "radio mode" and it counteracts the cooling flows of cold gas in undisturbed dark-matter haloes hosting galaxy clusters, which would otherwise show much higher star-formation of their central object than is observed. Previous work by Croton 2006 and De Lucia 2006 has shown that with the new semi-analytic model the population of local galaxies can be reproduced quite accurately. In order to study the evolution of the population out to higher redshifts, the semi-analytic predictions have been compared to a number of observations in various filter bands, in particular to two recent efforts to get a comprehensive multi-wavelength dataset of high redshift galaxies carried out by the DEEP2 (Davis 2001) and COSMOS (Scoville 2006) collaborations. The approach taken was to perform as broad a comparison as possible to gain firm constraints on the assumed physics in our model. Therefore a multitude of observational properties was contrasted with the model predictions such as clustering, luminosity functions, stellar mass functions, number counts per area and redshift to a certain magnitude limit. In order to facilitate the comparison between simulations and recent intermediate and high-redshift surveys, it is very useful to have a number of independent mock observations of the simulated galaxies, which provide good enough statistics to get a handle on cosmic variance. To this end I have devised a computer program that calculates the simulated galaxies lying on the backward light cone of a hypothetical observer out to arbitrarily high redshifts, taking advantage of the periodicity of the simulation box but avoiding replications. The output provides accurately interpolated redshifts, positions, observer frame and rest-frame magnitudes, dust extinction, as well as all the intrinsic galaxy properties like stellar mass and star formation rate. Utilising this tool it is also possible to make predictions for future galaxy surveys, deeper in magnitude and redshift than current ones. Presently the mock catalogues are used by the DEEP2 and COSMOS teams as a comparison sample in general and as a means to assess their selection effects and improve their data reduction in particular. First comparisons of counts in apparent magnitude and redshift gave promising results, showing good agreement in the low and intermediate range. The same holds for the angular clustering analysis except for the faintest magnitudes. Thus we conclude that our current understanding of the processes governing galaxy formation and evolution from the very first objects to the present day population is realistic but still incomplete. In particular the treatment of the interplay between star formation and negative feedback and the various processes influencing satellite galaxies in big galaxy clusters have potential for improvement. In the following I will give a brief outline of the thesis. After setting the stage for any kind of model in Chapter 1 by defining the geometry of the universe and the cosmological parameters that determine it, I will describe our semi-analytical model of galaxy formation in Chapter 2, where it will be also explained how to construct realistic mock observations of the simulated galaxies. First in Chapter 3 it will be verified that a simple model which assumes that galaxies are conserved but evolve in luminosity due to their star formation histories cannot account for the observed evolution of the galaxy population in the universe. This fact can be understood in the context of hierarchical models where massive and luminous galaxies assembled from smaller objects. Chapter 4 proceeds with exploring the predictions from the considerably more sophisticated semi-analytic model based on an N-body simulation of the hierarchical growth of dark matter structures. For this analysis a set of mock light-cones was constructed for direct comparison with the data which shows reasonably good agreement between model and observations at low redshift and for bright apparent magnitudes. These light-cones represent one of the largest samples of realistic mock observations currently available. They can be used for testing data analysis techniques usually applied to real observations on a well defined sample of artificial galaxies to verify how well the derivation of galaxy properties from the data works. In Chapter 5 we will demonstrate how one can measure the evolution of the galaxy merger rate from observing close projected galaxy pairs. Interestingly we find that the calibration needed for the conversion is significantly different from what has typically been assumed in previous studies. Additionally we will demonstrate that galaxy merger rates and dark-matter merger rates show considerably different evolution with redshift. Consequently we conclude that merger rate studies are less suitable as a probe of cosmic structure formation than initially assumed, but nonetheless they can be of great help to understand the formation and evolution of galaxies in a hierarchical universe. Finally these results will be summarised and discussed in Chapter 6 where I will also give a brief outlook on the future of this work, a short glimpse of which is already presented in the Appendix

Emergence of rich-club topology and coordinated dynamics in development of hippocampal functional networks in vitro.

Recent studies demonstrated that the anatomical network of the human brain shows a "rich-club" organization. This complex topological feature implies that highly connected regions, hubs of the large-scale brain network, are more densely interconnected with each other than expected by chance. Rich-club nodes were traversed by a majority of short paths between peripheral regions, underlining their potential importance for efficient global exchange of information between functionally specialized areas of the brain. Network hubs have also been described at the microscale of brain connectivity (so-called "hub neurons"). Their role in shaping synchronous dynamics and forming microcircuit wiring during development, however, is not yet fully understood. The present study aimed to investigate the role of hubs during network development, using multi-electrode arrays and functional connectivity analysis during spontaneous multi-unit activity (MUA) of dissociated primary mouse hippocampal neurons. Over the first 4 weeks in vitro, functional connectivity significantly increased in strength, density, and size, with mature networks demonstrating a robust modular and small-world topology. As expected by a "rich-get-richer" growth rule of network evolution, MUA graphs were found to form rich-clubs at an early stage in development (14 DIV). Later on, rich-club nodes were a consistent topological feature of MUA graphs, demonstrating high nodal strength, efficiency, and centrality. Rich-club nodes were also found to be crucial for MUA dynamics. They often served as broker of spontaneous activity flow, confirming that hub nodes and rich-clubs may play an important role in coordinating functional dynamics at the microcircuit level.M.S.S. is supported by a PhD studentship funded by a Core Award from the Medical Research Council and the Wellcome Trust to the Behavioural and Clinical Neuroscience Institute (MRC Ref G1000183; WT Ref 093875/Z/10/Z) and by the Studienstiftung des deutschen Volkes. Additional support for this study from the Biotechnology and Biological Sciences Research Council (BBSRC Ref BB/H008608/1) is gratefully acknowledged.This is the final published version. It first appeared at http://www.jneurosci.org/content/35/14/5459.full

Encoding cortical dynamics in sparse features.

Distributed cortical solutions of magnetoencephalography (MEG) and electroencephalography (EEG) exhibit complex spatial and temporal dynamics. The extraction of patterns of interest and dynamic features from these cortical signals has so far relied on the expertise of investigators. There is a definite need in both clinical and neuroscience research for a method that will extract critical features from high-dimensional neuroimaging data in an automatic fashion. We have previously demonstrated the use of optical flow techniques for evaluating the kinematic properties of motion field projected on non-flat manifolds like in a cortical surface. We have further extended this framework to automatically detect features in the optical flow vector field by using the modified and extended 2-Riemannian Helmholtz-Hodge decomposition (HHD). Here, we applied these mathematical models on simulation and MEG data recorded from a healthy individual during a somatosensory experiment and an epilepsy pediatric patient during sleep. We tested whether our technique can automatically extract salient dynamical features of cortical activity. Simulation results indicated that we can precisely reproduce the simulated cortical dynamics with HHD; encode them in sparse features and represent the propagation of brain activity between distinct cortical areas. Using HHD, we decoded the somatosensory N20 component into two HHD features and represented the dynamics of brain activity as a traveling source between two primary somatosensory regions. In the epilepsy patient, we displayed the propagation of the epileptic activity around the margins of a brain lesion. Our findings indicate that HHD measures computed from cortical dynamics can: (i) quantitatively access the cortical dynamics in both healthy and disease brain in terms of sparse features and dynamic brain activity propagation between distinct cortical areas, and (ii) facilitate a reproducible, automated analysis of experimental and clinical MEG/EEG source imaging data

Maturation trajectories of cortical resting-state networks depend on the mediating frequency band

The functional significance of resting state networks and their abnormal manifestations in psychiatric disorders are firmly established, as is the importance of the cortical rhythms in mediating these networks. Resting state networks are known to undergo substantial reorganization from childhood to adulthood, but whether distinct cortical rhythms, which are generated by separable neural mechanisms and are often manifested abnormally in psychiatric conditions, mediate maturation differentially, remains unknown. Using magnetoencephalography (MEG) to map frequency band specific maturation of resting state networks from age 7 to 29 in 162 participants (31 independent), we found significant changes with age in networks mediated by the beta (13–30 Hz) and gamma (31–80 Hz) bands. More specifically, gamma band mediated networks followed an expected asymptotic trajectory, but beta band mediated networks followed a linear trajectory. Network integration increased with age in gamma band mediated networks, while local segregation increased with age in beta band mediated networks. Spatially, the hubs that changed in importance with age in the beta band mediated networks had relatively little overlap with those that showed the greatest changes in the gamma band mediated networks. These findings are relevant for our understanding of the neural mechanisms of cortical maturation, in both typical and atypical development.This work was supported by grants from the Nancy Lurie Marks Family Foundation (TK, SK, MGK), Autism Speaks (TK), The Simons Foundation (SFARI 239395, TK), The National Institute of Child Health and Development (R01HD073254, TK), National Institute for Biomedical Imaging and Bioengineering (P41EB015896, 5R01EB009048, MSH), and the Cognitive Rhythms Collaborative: A Discovery Network (NFS 1042134, MSH). (Nancy Lurie Marks Family Foundation; Autism Speaks; SFARI 239395 - Simons Foundation; R01HD073254 - National Institute of Child Health and Development; P41EB015896 - National Institute for Biomedical Imaging and Bioengineering; 5R01EB009048 - National Institute for Biomedical Imaging and Bioengineering; NFS 1042134 - Cognitive Rhythms Collaborative: A Discovery Network

Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks.

Human functional magnetic resonance imaging (fMRI) brain networks have a complex topology comprising integrative components, e.g. long-distance inter-modular edges, that are theoretically associated with higher biological cost. Here, we estimated intra-modular degree, inter-modular degree and connection distance for each of 285 cortical nodes in multi-echo fMRI data from 38 healthy adults. We used the multivariate technique of partial least squares (PLS) to reduce the dimensionality of the relationships between these three nodal network parameters and prior microarray data on regional expression of 20 737 genes. The first PLS component defined a transcriptional profile associated with high intra-modular degree and short connection distance, whereas the second PLS component was associated with high inter-modular degree and long connection distance. Nodes in superior and lateral cortex with high inter-modular degree and long connection distance had local transcriptional profiles enriched for oxidative metabolism and mitochondria, and for genes specific to supragranular layers of human cortex. In contrast, primary and secondary sensory cortical nodes in posterior cortex with high intra-modular degree and short connection distance had transcriptional profiles enriched for RNA translation and nuclear components. We conclude that, as predicted, topologically integrative hubs, mediating long-distance connections between modules, are more costly in terms of mitochondrial glucose metabolism.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.PEV is supported by an MRC Bioinformatics Research Fellowship (MR/K020706/1). Functional MRI data acquisition was supported by a strategic award from the Wellcome Trust to the University of Cambridge (IMG, PBJ, ETB) and University College London (RJD, PF): the Neuroscience in Psychiatry Network (NSPN). Additional support was provided by the NIHR Cambridge Biomedical Research Centre. Access to gene expression data was provided by the Allen Institute for Brain Sciences Website: © 2015 Allen Institute for Brain Science. Allen Human Brain Atlas [Internet]. Available from: http://human.brain-map.org. FV is supported by a Gates Cambridge PhD studentship. KW is supported by the University of Cambridge MB/PhD Programme and the Wellcome Trust. We thank Gita Prabu, Roger Tait, Cinly Ooi, John Suckling and Becky Inkster for fMRI data collection and storage. ETB is employed half-time by the University of Cambridge and half-time by GlaxoSmithKline(GSK).This is the final version of the article. It first appeared from Royal Society Publishing via http://dx.doi.org/10.1098/rstb.2015.036

Specific Frontostriatal Circuits for Impaired Cognitive Flexibility and Goal-Directed Planning in Obsessive-Compulsive Disorder: Evidence From Resting-State Functional Connectivity.

BACKGROUND: A recent hypothesis has suggested that core deficits in goal-directed behavior in obsessive-compulsive disorder (OCD) are caused by impaired frontostriatal function. We tested this hypothesis in OCD patients and control subjects by relating measures of goal-directed planning and cognitive flexibility to underlying resting-state functional connectivity. METHODS: Multiecho resting-state acquisition, combined with micromovement correction by blood oxygen level-dependent sensitive independent component analysis, was used to obtain in vivo measures of functional connectivity in 44 OCD patients and 43 healthy comparison subjects. We measured cognitive flexibility (attentional set-shifting) and goal-directed performance (planning of sequential response sequences) by means of well-validated, standardized behavioral cognitive paradigms. Functional connectivity strength of striatal seed regions was related to cognitive flexibility and goal-directed performance. To gain insights into fundamental network alterations, graph theoretical models of brain networks were derived. RESULTS: Reduced functional connectivity between the caudate and the ventrolateral prefrontal cortex was selectively associated with reduced cognitive flexibility. In contrast, goal-directed performance was selectively related to reduced functional connectivity between the putamen and the dorsolateral prefrontal cortex in OCD patients, as well as to symptom severity. Whole-brain data-driven graph theoretical analysis disclosed that striatal regions constitute a cohesive module of the community structure of the functional connectome in OCD patients as nodes within the basal ganglia and cerebellum were more strongly connected to one another than in healthy control subjects. CONCLUSIONS: These data extend major neuropsychological models of OCD by providing a direct link between intrinsically abnormal functional connectivity within dissociable frontostriatal circuits and those cognitive processes underlying OCD symptoms.This research was funded by a Wellcome Trust Senior Investigator Award (104631/Z/14/Z) awarded to T.W. Robbins. Work was completed at the Behavioural and Clinical Neuroscience Institute, University of Cambridge, Cambridge, UK, supported by a joint award from the Medical Research Council and Wellcome Trust (G00001354). M.M. Vaghi is supported by a Pinsent Darwin Scholarship in Mental Pathology and a Cambridge Home and EU Scholarship Scheme (CHESS) studentship. P.E. Vértes is supported by the Medical Research Council (grant no. MR/K020706/1). A.M. Apergis-Schoute is supported by the Wellcome Trust above. V. Voon is a Wellcome Trust Fellow

Somatosensory cortex functional connectivity abnormalities in autism show opposite trends, depending on direction and spatial scale

Functional connectivity is abnormal in autism, but the nature of these abnormalities remains elusive. Different studies, mostly using functional magnetic resonance imaging, have found increased, decreased, or even mixed pattern functional connectivity abnormalities in autism, but no unifying framework has emerged to date. We measured functional connectivity in individuals with autism and in controls using magnetoencephalography, which allowed us to resolve both the directionality (feedforward versus feedback) and spatial scale (local or long-range) of functional connectivity. Specifically, we measured the cortical response and functional connectivity during a passive 25-Hz vibrotactile stimulation in the somatosensory cortex of 20 typically developing individuals and 15 individuals with autism, all males and right-handed, aged 8-18, and the mu-rhythm during resting state in a subset of these participants (12 per group, same age range). Two major significant group differences emerged in the response to the vibrotactile stimulus. First, the 50-Hz phase locking component of the cortical response, generated locally in the primary (S1) and secondary (S2) somatosensory cortex, was reduced in the autism group (P < 0.003, corrected). Second, feedforward functional connectivity between S1 and S2 was increased in the autism group (P < 0.004, corrected). During resting state, there was no group difference in the mu-α rhythm. In contrast, the mu-β rhythm, which has been associated with feedback connectivity, was significantly reduced in the autism group (P < 0.04, corrected). Furthermore, the strength of the mu-β was correlated to the relative strength of 50 Hz component of the response to the vibrotactile stimulus (r = 0.78, P < 0.00005), indicating a shared aetiology for these seemingly unrelated abnormalities. These magnetoencephalography-derived measures were correlated with two different behavioural sensory processing scores (P < 0.01 and P < 0.02 for the autism group, P < 0.01 and P < 0.0001 for the typical group), with autism severity (P < 0.03), and with diagnosis (89% accuracy). A biophysically realistic computational model using data driven feedforward and feedback parameters replicated the magnetoencephalography data faithfully. The direct observation of both abnormally increased and abnormally decreased functional connectivity in autism occurring simultaneously in different functional connectivity streams, offers a potential unifying framework for the unexplained discrepancies in current findings. Given that cortical feedback, whether local or long-range, is intrinsically non-linear, while cortical feedforward is generally linear relative to the stimulus, the present results suggest decreased non-linearity alongside an increased veridical component of the cortical response in autism

Conservative and disruptive modes of adolescent change in human brain functional connectivity.

Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: "conservative" and "disruptive." Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman's correlation between edgewise baseline FC (at 14 y, [Formula: see text]) and adolescent change in FC ([Formula: see text]), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas

Two human brain systems micro-structurally associated with obesity

The relationship between obesity and brain structure is incompletely understood. Using diffusion-weighted MRI from ∼30,000 UK Biobank participants we test the hypothesis that obesity (waist-to-hip ratio, WHR) is associated with regional differences in two micro-structural MRI metrics: isotropic volume fraction (ISOVF), an index of free water, and intra-cellular volume fraction (ICVF), an index of neurite density. We observed significant associations with obesity in two coupled but distinct brain systems: a prefrontal-temporal-striatal system associated with ISOVF and a medial temporal-occipital-striatal system associated with ICVF. The ISOVF∼WHR system colocated with expression of genes enriched for innate immune functions, decreased glial density, and high mu opioid (MOR) and other neurotransmitter receptor density. Conversely, the ICVF∼WHR system co-located with expression of genes enriched for G-protein coupled receptors and decreased density of MOR and other receptors. To test whether these distinct brain phenotypes might differ in terms of their underlying shared genetics or relationship to maps of the inflammatory marker C-reactive Protein (CRP), we estimated the genetic correlations between WHR and ISOVF (rg = 0.026, P = 0.36) and ICVF (rg = 0.112, P < 9 × 10−4) as well as comparing correlations between WHR maps and equivalent CRP maps for ISOVF and ICVF (p<0.05). These correlational results are consistent with a two-way mechanistic model whereby genetically determined differences in neurite density in the medial temporal system may contribute to obesity, whereas water content in the prefrontal system could reflect a consequence of obesity mediated by innate immune system activation

Maturation trajectories of cortical resting-state networks depend on the mediating frequency band.

The functional significance of resting state networks and their abnormal manifestations in psychiatric disorders are firmly established, as is the importance of the cortical rhythms in mediating these networks. Resting state networks are known to undergo substantial reorganization from childhood to adulthood, but whether distinct cortical rhythms, which are generated by separable neural mechanisms and are often manifested abnormally in psychiatric conditions, mediate maturation differentially, remains unknown. Using magnetoencephalography (MEG) to map frequency band specific maturation of resting state networks from age 7 to 29 in 162 participants (31 independent), we found significant changes with age in networks mediated by the beta (13-30 Hz) and gamma (31-80 Hz) bands. More specifically, gamma band mediated networks followed an expected asymptotic trajectory, but beta band mediated networks followed a linear trajectory. Network integration increased with age in gamma band mediated networks, while local segregation increased with age in beta band mediated networks. Spatially, the hubs that changed in importance with age in the beta band mediated networks had relatively little overlap with those that showed the greatest changes in the gamma band mediated networks. These findings are relevant for our understanding of the neural mechanisms of cortical maturation, in both typical and atypical development